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Recent experimental and theoretical studies show that energy efEciency which measures the amount of infor¬ 
mation processed by a neuron with per unit of energy consumption, plays an important role in the evolution of 
neural systems. Here, we calculated the information rates and energy efficiencies of the Hodgkin-Huxley (HH) 
neuron model at different temperatures in a noisy environment. We found that both the information rate and 
energy efficiency are maximized by certain temperatures. Though the information rate and energy efficiency can¬ 
not be maximized simultaneously, the neuron holds a high information processing capacity at the temperature 
corresponding to maximal energy efEciency. Our results support the idea that the energy efEciency is a selective 
pressure that inEuences the evolution of nervous systems. 
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Information processing in the nervous system is 
metabolically expensive. Human brain is only about 
2% of the total body weight, but consumes about 20% 
of the resting metabolic energy. [^1 The large metabolic 
energy requirement of nervous system could constrains 
the size and structure of the brain, and may have 
largely optimized the nervous system by favouring en¬ 
ergy efficient neural morphologies, codes, wiring pat¬ 
terns and brain structures.The energy efficiency of 
nervous system have possibly been greatly optimized 
by natural selection. 

There are many factors that influence the energy 
efficiency. Of all the information processing activities, 
action potential (AP) makes a great contribution of 
total consumed energy. AP on its own does not 
require energy, however, restoring the transmembrane 
ionic concentration gradients through the Na^/K'^ 
ion pump is an energy consuming process which re¬ 
lies on the energy released by adenosine triphosphate 
(ATP) hydrolysis. In the classic Hogkin-Huxley (HH) 
model, due to overlap between Na^ entry and 
outflow at the same time, the flow of Na^ and 
ions during AP largely exceed the minimum required 
by a pure capacitor, and thus waste large amounts 
of energy to restore the membrane potential.Re¬ 
cent investigations on the nonmyelinated mossy fibers 
of the rat hippocampus shows that minimizing the 
overlap of these two ion fluxes improves the energy ef¬ 
ficiency of AP generation. Energy efficiencies of the 
neural systems are also constrained by their size. Re¬ 


cent theoretical and numerical studies suggest that in 
the noisy environment the energy efficiency is max¬ 
imized by the number of ion channels in a single 
neuron, or the number of neurons in a neuronal 
circuits. 

Temperature is an important factor that constrain 
the energy efficiency of neural systems. It influences 
the conductance, activation and inactivation of all 
ion channels as a global perturbation. Tempera¬ 
ture sensitivities of ion channels present a challenge 
to maintaining stable functions over an extended tem¬ 
perature range.Yu et ai, found that warmer body 
temperatures reduce the Na ^/overlap and facili¬ 
tate single energy efficient action potentials. 

However, information in neural system is often rep¬ 
resented in groups of APs generated by neuron popu¬ 
lation or a single neuron within a time period, rather 
than a single AP. Therefore, it is interesting to study 
the effects of temperature on the energy efficiency 
in processing input signals with the consideration of 
noise perturbation. In this research, we stimulate sin¬ 
gle compartment HH neuron with synaptic inputs in 
presence of Gaussian white noise. The effective energy 
efficiency rate was introduce to measure how much in¬ 
formation is encoded in a unit time using unit energy 
in response to the synaptic stimuli. We then investi¬ 
gated the influence of temperature on the energy ef¬ 
ficiency of information processing by the HH neuron 
and found that the energy efficiencies exhibit a maxi¬ 
mum at certain temperatures. 
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Figure 1: The information capacity of the HH neuron as a function of temperature for different input signal 
strength (unit: /iA/cm^). (a) total entropy rate, (b) noise entropy rate, and (c) the information rate. 


The single compartment HH neuron with Gaussian 
white noise ^(t) and external current injection Ie(t) is 
described by following equations. 

= —gNa'm^h{V — Vno) — gKn^{V — Vk) 
at 

-gL{V-VL)+m + Ie{t), ( 1 ) 

dx 

— = aa;{y){l - x) - PxiV) {x = m,h,n).{2) 

at 


tion/inactivation variables m, h or n. The voltage- 
dependent rate functions of gating states transition 
ax{V) and l3a:{V) are 

am = = ((i(T)4.0e“^; 

an = /3„ = (;!)(r)0.125e“^; 

tth = (l)iT)0.07e~^, l3h = <l}iT)^^ -. 

e 10 -1-1 


Eq. (1) shows the varations in membrane potential (H) 
over time where Cm = l/rF/cm^ is the membrane 
capacitance. The maximum channel conductances 
and reversal potentials are typical values: gNa = 
120mS/cm^, gx = 36.0mS/cm^, = 0.33mS/cm^,k 

VNa = 50mV, Vk = -77mV and Vl = -54.4mV. 
m^h and are ratios of sodium channels and potas¬ 
sium channels at open state to all sodium channels 
and potassium channels, respectively, ^{t) is the zero- 
mean Gaussian white noise with 


< m >= 0; < mat') >= 2D5{t - t'), (3) 


where D is the noise intensity. Ie{t) = '^jliynit) is 
composed of many single synaptic-current-like pulses 
liyn{t) , a single pulse have the form 


Ilynit) = 


Io{t - ti)e ti<t<tl, 

0, t <ti or t > ti, 


(4) 


where j denote the jth pulse, Iq represents the 
strength of the synaptic-current-like pulse, t^, t^ = 
tl + At [At = Sms in our simulation) are the time 
when released neural transimitters first contact with 
the post-synaptic membrane (the pulse start) and the 
cut-off time of the pulse (the pulse end), represen¬ 
tatively. T = 2ms is the time delay constant. The 
spiking time of pulses follow Poisson distribution with 
the average time interval 100ms as expected for a real 
neuron. 

Eq. (2) represents three similar gating equations 
of Na'^ and K~^ channels, x can be any activa¬ 


where <?!>(T) is the temperature dependent function^^^l 
(/>(T) = (5) 

The equations are numerically integrated using 
Euler method with time step dt = O.Olms. The total 
simulation time lasted for 3600s in order to estimate 
the information entropy (see below). 

To measure the average information processing ef¬ 
ficiency of our HH neuron under the mimicked ex¬ 
ternal synaptic pulses, we calculated the information 
entropy rate and information rate using the method 
of Strong et al^ 

We first calculated the entropy rate of the spike 
train. First, we divided the spike train in response to 
set stimulus into bins of size At with value 1 (AP) or 
0 (no AP). Then we divided the bin train into T / At 
words. Gonsidering the size of data set, the total en¬ 
tropy rate is defined as 

Htotai = -T lim ^'^Pi\og 2 Pi, (6) 

i 

where Pi represents the probability of *-th word. Sim¬ 
ilar to the total entropy rate, we evaluate the infor¬ 
mation rate interrupted by noise using average noise 
entropy rate 

Hno%se=<- lim Pi{t)\og 2 Pi{t) >t, (7) 

T-s-oo 1 ^' 
i 

where Pi{t) is the time-dependent word probability. 
We then calculated information rate, which measures 
the average information encoded into the spike trains 
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Figure 2: Energy consumption rates as a function of temperatures of the HH neuron, (a) total energy consump¬ 
tion rates; (b) energy consumption rates induced by noise alone; (c) the energy cost of a single AP. The legend 
of A and C denote the stimus strength with unit /iA/cm^. 


in response to stimulus, by subtracting noise entropy 
rate from total entropy rate. 

I — titotal Hnoise- ( 8 ) 

Fig. 1 (a) , (b) and (c) demonstrate the total en¬ 
tropy rate, noise entropy rate, and information rate as 
functions of temperature for HH neuron in response 
to signals with different strength, respectively. We 
observed that both the total entropy rate and noise 
entropy rate decrease as the temperature increases. 
However, when temperature is high, the total entropy 
rate decreases rapidly, whereas the noise entropy rate 
shows a local maxima. The information rate, i.e. the 
difference between total entropy rate and noise en¬ 
tropy rate, exhibits a global maxima at certain tem¬ 
peratures between 0 to 25°C. 

As shown in Fig. 1, for strong signals, the total en¬ 
tropy rate is high but the noise entropy is low. This 
may be because strong signals can induce more APs, 
which in turn increases the total information received 
by the neuron (Fig. 1 (a)). Strong signals also have 
high probability to induce APs, which increases the 
coherence of firings across the trials, thus decrease the 
noise entropy in the spike trains (Fig. 1 (b)). As a re¬ 
sult, the information rate is high for strong signals and 
decreases as signal strength decrease (Fig. 1 (c)). 

We also calculated the energy consumption rate of 
the HH neuron in the above signal processing event. 
The energy cost of a neuron are calculated mostly by 
integrating the Na'^ or K'^ current over time, and 
then convert them into the number of Na^ or 
ions. Since the Na^/K'^ pump hydrolyzes one ATP 
molecule for three Na^ ions extruded and two K'^ 
ions imported, the energy cost is estimated by the 
amount of ATP molecules are expended, or in units of 
energy if the powers of the Na ^/pump are mea¬ 
sured in real experiment (for example, 50kJ/mol in the 
heart Recently, Moujahid et al. proposed a novel 
method to estimate the energy cost directly from the 


equivalent electrical circuit of the HH neuron. 

Here in this paper we use a similar but more direct 
method to calculate the energy cost of the above HH 
neuron in processing signals. We consider the energy 
provided by or leak out to external environment. In 
the equivalent electrical circuit of HH neuron, the en¬ 
ergy is consumed by three types of ion channel con¬ 
ductance and can be provided or consumed by exter¬ 
nal electrode-clamping circuits. Thus the total energy 
consumption rate is equal to the sum of all four powers 

= Pe{t) +PNai.t) +PK{t) +pL{t) (9) 

= VIe{t)+iNaiV -VNa) 

+iK{V-VK)+iL{V-VL) ( 10 ) 

= VIeit) - gNam^HV -VNaf 
-gKn\V - Vk? - gdv - 

In Eq. (11), the minus sign in front of ion chan¬ 
nel terms represent energy consumption while ignored 
plus sign in front of external stimuli term represents 
power supply. To calculate total energy consumption, 
we reversed the signs on the right hand side of Eq. (11) 
to denote energy consumption of ion channels as posi¬ 
tive and power supplied (or consumption if the value of 
VIe{t) itself is negative) by external circuits as neg¬ 
ative. To estimate Etotai, the total energy cost per 
unit time in the signal detection progress, dE{t)/dt is 
integrated over time and then divided by total sim¬ 
ulation time. The energy cost rate induced by noise 
alone (Enoise) is calculated as above but without in¬ 
put signals. The energy cost of a single AP {Esingle) 
is calculated by integrating dE{t)/dt in the period of 
the AP. 

As demonstrated in Fig. 2 (a) , the total energy 
consumption rate decreases as temperature increases. 
At lower temperature, the energy consumption rate 
is almost identical for different signal strengths. At 
high temperature, the neuron consumes more energy 
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to processing strong signals. The dependence of Enoise 
on temperature decreases with temperature and has 
two distinct stages: at lower temperature (less than 
15°C) it decreases almost lineally and at higher tem¬ 
perature (above 15°C) it almost flattens (Fig. 2 (b)). 
As displayed in Fig. 2(c), the average energy cost of 
single AP is identical to different signal strength, be¬ 
cause the shape of AP almost does not change as sig¬ 
nal changes. As the temperature increases, the energy 
cost of an AP decreases, which is in line with previous 
studies. 



Figure 3: Energy efficient as a function of temperature 
for different input signal strengths. The solid markers 
denote the maximum energy efficiency. 

Finally, we define the energy efficiency as the in¬ 
formation rate divided by energy consumption rate, 

i.e. 



which measures how much information is effectively 
processed by the neuron in response to the stimulus 
and consuming one unit of energy. As shown in Fig. 3, 
with the increase in temperature, the energy efficiency 
first increases and then decreases, passing through a 
global maximum as a function of temperature for dif¬ 
ferent signal strength. The existence of maximum im¬ 
plies that the temperature could maximize the energy 
efficiency of the neuron in the processing of input sig¬ 
nals. Meanwhile, with the increasing signal strength, 
the energy efficiency increases, as well as the temper¬ 
ature corresponding to maximal energy efficiency. 

Comparing data at the optimal temperature, i.e. 
the temperature with maximum energy efficiency, and 
at the temperature with maximum information rate 
(see Table. 1, variables at the optimal temperature are 
marked with superscript ‘o’), we found that the tem¬ 
peratures that maximize energy efficiency {T°) and 


that which maximize information rate (T^) for a par¬ 
ticular signal strength is not the same; information 
rate and energy efficiency cannot be maximized by the 
temperature simultaneously. However, at the temper¬ 
ature the energy efficiency is maximized, the corre¬ 
sponding information rate is in the range of 70% to 
90% of the maximal information rate (j4^). The 
signal detection rate PdetecUon^ which is ratio of the 
number of APs induced immediately (less thant Sms) 
after signals are applied and the total number of ap¬ 
plied signals at T°, is in the range from 80% to 99%. 
Thus, it is concluded that at the temperature the en¬ 
ergy efficiency is maximized, the neuron can still keep 
relatively high capacity to process information. At the 
same time, at the optimal temperature T°, the spon¬ 
taneous firing rate is very low < 0.27Hz ), and 

the noise induced energy cost is less than 40% of total 
energy cost ). 

Table 1. Variables at optimal temperature 


I(^) r°(°c) r^(°c) 


po 

detection 


/noise (Hz) 


^total 


6 

14 

6 

0.7094 

0.7922 

0.262 

0.3012 

8 

16 

10 

0.8341 

0.9367 

0.075 

0.2943 

10 

18 

12 

0.8509 

0.9766 

0.018 

0.3219 

12 

20 

12 

0.8998 

0.9894 

0.005 

0.3573 

14 

22 

14 

0.8898 

0.9941 

0 

0.3954 


In conclusion, we investigated the influence of tem¬ 
perature on the capacity of information processing and 
energy efficiency of the classic HH neuron in response 
to synaptic inputs. We found that temperature can 
maximize the information capacity, as well as the en¬ 
ergy efficiency of HH neuoron. The difference between 
the temperature that maximizes energy efficiency and 
that which maximizes information rate for a partic¬ 
ular signal strength produces a conflict; information 
rate and energy efficiency cannot be maximized simul¬ 
taneously. This conflict has been observed in models 
of spiking neurons and neural codes,adding 
to the numerous lines of evidence suggesting that en¬ 
ergy is a selective pressure that has influenced the evo¬ 
lution of neural systems.Our study further show 
that at the temperature the energy efficiency is max¬ 
imized, the neuron can still keep its information pro¬ 
cessing capacity high, comparing its values at maxi¬ 
mum. Considering that the typical living temperature 
of Loligo (the squid used by Hogkin and Huxleyl^"^!) is 
10-26‘^C[22.23]^ 

our result of optimized temperature is 
consistant with real data. 
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